Theta-Induced Electric Dipole Moment of the Neutron via QCD Sum Rules 
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Using the QCD sum rule approach, we calculate the electric dipole moment of the neutron induced 
by a vacuum S-angle to approximately 40-50% precision, dn = 2.4 x 10~^^6e ■ cm. Combined with 
the new experimental bound, this translates into the limit |^| < 3 x 10"^". 
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The impressive experimental limits on the electric 
dipole moments (EDMs) of neutrons and heavy atoms in 
general put a very strong constraint on possible flavor- 
conserving CP-violation around the electroweak scale [1] . 
This precision means that EDMs can in principle probe a 
high energy scale by limiting the coefficients of operators 
with dimension> 4, such as GG, GGG, and qGaj^q etc.. 
However, in practice, while these operators can be per- 
turbatively evolved down to a scale of order 1 GeV, the 
ultimate connection between high energy parameters and 
low energy EDM observables necessarily involves non- 
perturbative physics. 

The 9-teim, OGG, which has dimension=4 may be 
interpreted as the lowest dimension CP-violating oper- 
ator, un-suppressed by the high energy scale. Exper- 
imental tests of CP symmetry suggest that 6 is small 
and, among different CP-violating observables, the EDM 
of the neutron (c?„(0)) is the most sensitive to its value 
[3, 4]. However, the calculation of dn{0) is a long stand- 
ing problem [5-7] . According to Ref. [6] , an estimate can 
be obtained within chiral perturbation theory, relying on 
the numerical dominance of contributions proportional 
to lnm,r near the chiral limit. However, there are also 
non-logaritmic contributions, incalculable within this for- 
malism, which in principle can be numerically more im- 
portant than the logarithmic piece away from the chiral 
limit. As a consequence one is unable to estimate the 
uncertainty of the prediction [6] . 

There are a number of important incentives for refin- 
ing the calculation of d„ (6) . While 6 primarily arises as a 
fundamental vacuum parameter, one may also induce cal- 
culable corrections via the integration over heavy fields, 
and such corrections are significant in theories where the 
fundamental parameter is absent. One scenario of this 
type arises as a result of exact P or CP symmetries [8], 
whose subsequent spontaneous breakdown allows radia- 
tive corrections to induce an effective 0-parameter, which 
is then the dominant source for ci„. Alternatively, the 
fundamental 6 parameter may be removed via the axion 
mechanism, but then a small but finite effective 6'-angle 
will again survive, induced by the linear shift of the ax- 
ionic potential due to higher-dimensional CP-odd opera- 
tors [9] . In this case, the calculation of dn (0) is part of a 
more generic calculation of dn- 

In this letter, we apply the QCD sum rule method [2] to 



obtain an estimate for d„(0) beyond chiral perturbation 
theory. Within currently available analytic techniques, 
QCD sum rules appears the most promising approach to 
this problem as it has, in particular, been used success- 
fully in the calculation of certain baryonic electromag- 
netic form factors [10, 11]. Within the sum rule formal- 
ism, physical properties of the hadronic resonances are 
expressed via a combination of perturbative and nonper- 
turbative contributions, the latter parametrized in terms 
of vacuum quark and gluon condensates. We note that 
previously QCD sum rules were used to estimate the neu- 
tron EDM induced by a CP-odd color electric dipole mo- 
ment of quarks [12, 13] . Surprisingly, the results dis- 
agree with the naive estimates based on the chiral loop 
approach. The calculation of dn{0) using QCD sum rules 
will certainly help to resolve this controversy. 

The approach we shall use follows recent work [14] on 
the 6'-induced p-meson EDM in reducing the operator 
product expansion to a set of vacuum condensates taken 
in an electromagnetic and topologically nontrivial back- 
ground. Expansion to first order in 9 results in the ap- 
pearance of matrix elements which can be calculated via 
the use of current algebra [15, 16]. In this approach the 
^-dependence arises with the correct quark mass depen- 
dence, and the relation to the U(l) problem becomes ex- 
plicit as dn{0) vanishes when the mass of the U(l) "Gold- 
stone boson" is set equal to the mass of pion. 

The starting point for the calculation is the correlator 
of currents ry„(x) with quantum numbers of the neutron 
in a background with nonzero 9 and an electromagnetic 
field F^^, 

H(Q2) = , / d4^e^?-(Olr{r;„(a;)rJ„(0)}lO)e,F, (1) 



where = —q^, with q the current momentum. 

Firstly, it is important to note that if CP-symmetry 
is broken by a generic quark-gluon CP-violating source 
(^-term in our case), the coupling between the physical 
state (neutron) described by a spinor v and the current 
rjn acquires an additional phase factor 
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(2) 



The existence of this unphysical phase a is apparent al- 
ready when one considers the sum rule for the neutron 
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FIG. 1: Various contributions to the CP-odd structure {Fcr, /j}. 
(a) is the leading order contribution while (b), (c) and (d) con- 
tribute at subleading order. 



mass, which in the absence of CP-invariance can have an 
additional Dirac structure proportional to 175. When we 
turn to electromagnetic form factors, this angle can mix 
electric (d) and magnetic (/i) dipole moment structures 
and complicate the extraction of d from the sum rule. 
Consequently, in general the neutron double pole contri- 
bution on the phenomenological side of the sum rule will 
be proportional to the following expression 

Ua{^ + mn){^J.F(T - dFa){4 + ■mn)Ua 
= rrinlfiFa - dPa, ^} 

+ ^(At±ad)0(i^) +^((i±aA*)C(^'), (3) 

in which 0{F) and 0{F) are operators depending on the 
electromagnetic field strength F and its dual F which 
contain an even number of 7-matrices. We see that only 
Lorentz structures with an odd number of 7-matrices 
are independent of a. In calculating rf„, it is then clear 
that we should study the operator {Fa,c^, as this is the 
unique choice with a coefficient independent of a. 

The interpolating current r/„ is conveniently paramet- 
rised in the form, 



Vn 



where the two contributions are given by 

ji = 2eabc{dlC'y5Ub)dc 
h = 2eabc{dlCub)'y5dc- 



(4) 

(5) 
(6) 



j2 vanishes in the nonrelativistic limit, and lattice simu- 
lations have shown that ji indeed provides the dominant 
projection onto the neutron (see e.g. [17]). From (2), we 
may define the coupling to the neutron state in the form 
(0|r;„|iV) = (Ai+/3A2)t;. 

Within the sum rules formalism, one has the impera- 
tive of suppressing the contribution of excited states and 



higher dimensional operators in the operator product ex- 
pansion (OPE), and thus its convenient to choose /3 to 
this end. loffe has shown that /? = — 1 is an apparently 
optimal choice for the mass sum rule. liowevcr, it is 
clear that this optimization may differ for different phys- 
ical observables. We shall therefore keep /3 arbitrary, and 
optimize once we have knowledge of the structure of the 
sum rule. 

We now proceed to study the OPE associated with 
(1). The relevant diagrams we need to consider arc 
shown in Fig. 1 ((a), (b) and (c)). Diagrams of the form 
(d), although suffering no loop factor suppression, are 
nonetheless suppressed due to combinatorial factors and 
the small numerical size of {{qq))^- 

In parametrizing 0, we shall take a general initial con- 
dition in which a chiral rotation has been used to generate 
a 75 -mass, so that 



■ Oqiri: 



g/n5«/ + eGgG«,G», + ...,(7) 
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in which we restrict to qf = u, d, and so = 
iTiumd/ (ruu + ma)- The physical parameter is of course 
6 = Oq + 9o, but we shall keep the general form (7) and 
calculate the OPE as a function of both phases. The in- 
dependence of the final answer of 9q — 0c will provide a 
nontrivial check on the consistency of our approach. We 
shall find that this requires the consideration of mixing 
with additional currents, a point we shall come to shortly. 

The vacuum structure is conveniently encoded in a 
generalised propagator expanded in the background field 
and the associated condensates. We work as usual with 
a constant background electromagnetic field, so that 
Af^{x) = — ii^p^(0)x'^, and use a fixed point gauge for 
the gluon potential, A^^fix) = -^Gl^^iO^fx" . The elec- 
tromagnetic field dependence is determined in terms of 
the magnetic susceptibilities x, « and ^, defined as [10]: 



{0\q(Tu,^q\0)F 
gmGlX)q\0)F 

2g{0\q^s{GlJ'')q\0)F 



KqF^AO\qq\0) 
i^qF^u{0\qq\0), 



while the ^-dependence is either explicit in the case of 
6q, or extracted via use of the anomalous Ward identity 
(see e.g. [16]) in the case of Oq- The dependence of 
the OPE then follows as a consequence of and 
disappears when J7(l)-symmetry is restored (m^ — > ttItt) 
[14]. 

Defining iS{q) = {0\qa{x)q^{0)\0)F,e, and ignoring a 
trivial (5-function over colour indices, the leading order 
propagator adapted to the CP-odd sector and appropri- 
ate for Fig. 1(a), takes the form. 



Slo{x) 



(1 - i0qj5)ab 



-^^^^«/3a;„(7/375)6a - ^{Fa{l+ieG7,))ba, (8) 
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where we have denoted Xq = Xq {q<i) , and we shall hence- 
forth follow [10] and assume that Xq = X^q etc., with 
flavour independent parameters x, k, ^. 

Diagrams (b) and (c) in Fig. 1 require, in addition, 
the leading order expansion in the background gluon and 
electromagnetic fields, and Snlo is consequently more 
involved. However, since the approach is analogous to 
that for the leading order terms, we shall defer full de- 
tails [18], and simply present the final result for the OPE 
structure arising from these diagrams. In momentum 
space we find 



n(g') = 



9m^ 
64^ 



{qq){F(T, ^ 



x(/3+l)'(4ed-e„)ln 



-4(/3-l)2ed l + -(2« + e) 



1 



1 _i 

-| ((4/?2 - 4/3 + 2)ea + + 2/3 + l)e„) ^ 



A2 



(9) 



where the first, second, and third lines correspond to con- 
tributions from diagrams (a), (b)+(c), and (c), in Fig. 1, 
and is an infrared cutoff. The diagrams (d) which 
contribute at sublcading order in are more problem- 
atic because they involve correlators not calculable in the 
chiral approach, e.g. 



n 



/ d''xe^''^0\T{{qqf{x),m.^qfij,qf}\0) 
J ^ 



(10) 



We estimate such contributions via saturation with the 
physical rj meson. The result turns out to be parametri- 
cally smaller than any term listed in Eq. (9). 

At different stages of the calculation we observe the ap- 
pearance of the unphysical phase 9g — Oq which is closely 
related with the non-invariance of the currents ji and 
j2 under a chiral transformation. This invariance is, in 
fact, restored when we consider the mixing of ji and j2 
with another set of currents, ii = 2eabc{d^Cub)dc and 
«2 = '^^abc{d^'y5Cuh)-f5dc. The mixing between these two 
sets is explicitly proportional to 9q — 9q. When properly 
diagonalized on the sum rule for ^, the linear combination 
of these two sets of currents provides additional contri- 
butions proportional to 9g — 9q, exactly cancelling the 
unphysical piece of n(Q^). Consequently, it is equivalent 
to take 6q = 9g as the most convenient basis when work- 
ing with ji and j2 where the unwanted phase, and the 
mixing between ji^2 and ii,2, are simply absent. This 
situation resembles the problem in obtaining the "cor- 
rect" quark mass behavior for the EDM of p, addressed 
in detail in [14]. 

Inspection of the expression (9) shows that the stan- 
dard choice of /? = — 1 appropriate for CP-even sum rules 
will not be useful here as it removes the leading order 
contribution. In general there are two motivated tech- 
niques for fixing the mixing parameter jS: (1) at a local 
extremum; or (2) to minimize the effects of the contin- 
uum and higher dimensional operators. We find in this 



case that extremizing in the parameter /3 also leads to 
the unappealing cancelation of the leading order contri- 
bution. Thus the most natural procedure appears to be 
to choose /3 in order to cancel the subleading infrared 
logarithm which is ambiguous as a result of the required 
infrared cutoff. This procedure actually mimics the orig- 
inal motivation for /3 = — 1 in the CP-even case. We 
then take /3 = 1, and it is this choice that we shall now 
contrast with the phenomenological side of the sum rule. 
It is worth to remarking, however, that use of the "lat- 
tice" current with /? = will also produce a numerically 
similar result. 

On the phenomenological side of the sum rule we have 
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(11) 



We retain here double and single pole contributions, the 
latter corresponding to transitions between the neutron 
and excited states, but the exponentially suppressed con- 
tinuum contribution will be ignored as it is not, in this 
case, sign definite. Here we use A = Ai -|- (3X2 and A is 
an effective constant parametrizing the single pole con- 
tributions. 

After a Borel transform of (9) and (11) and using /3 = 1 
as discussed earlier, we obtain the sum rule 



Xnindn + AM = ~-^^^m*{^q)e 

1 



Ax{Aed - Bu) 



2M2 



C(4ed + 8e„) 



• (12) 



The coupling A present in (12) may be obtained from 
the well known sum rules for the tensor structures 1 and 
^ in the CP even sector (see e.g. [17] for a recent review). 
We shall construct two sum rules in this way. 

• (a) Firstly, we extract a numerical value for A via 
a direct analysis of the CP even sum rules. This 
analysis has been discussed before and will not be 
reproduced here (see e.g. [17]). One uses /3 = — 1, 
and obtains (27r)''A ~ 1.05 ± 0.1. 

• (b) As an alternative, we extract A explicitly as a 
function of /3 from the CP-even sum rule for ^, and 
substitute the result into (12) choosing /3 = 1. 

The conventional approach which we shall adopt here 
is to assume that A is independent of M, and thus the 
left hand side of (12) is linear in M"^ provided that A is 
constant in the appropriate region for the Borel param- 
eter. The latter point in case (b) above may be checked 
explicitly. 

It is convenient to define an additional function z/(M^), 



!/(M2) = 



'111 in . 



dn + 



AM^ 



(13) 



which is then determined by the right hand side of (12). 
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FIG. 2: The neutron EDM function v{M'^ {CfN'^)) is plotted ac- 
cording to the sum rules (a) and (b). The dashed line shows the 
contribution from the leading order term only. 



The two sum rules described above for Va and are 
plotted in Fig. 2, where the effect of the higher dimen- 
sional terms in (12) proportional to ^ is also displayed. 
viM"^) is to be interpreted as a tangent to the curves in 
Fig. 2. For numerical calculation we make use of the fol- 
lowing parameter values: For the quark condensate, we 
take (0|gg|0) = -(0.225 GeV)^, while for the condensate 
susceptibilities, we have the values x = — 5.7±0.6 GeV~^ 
[19], and i = -0.74 ± 0.2 [13]. Note that x, which 
enters at 0(1/M^), since it is dimensionful, is numeri- 
cally significantly larger than ^. In extracting A in case 
(b) we also set a relatively large continuum threshold 
So = (2GeV)^ for consistency with the CP-odd sum rule 
in which this continuum is ignored. 

One observes that both sum rules have extrema consis- 
tent to ^ 50%, suggesting that our procedure for fixing 
the parameter /? is appropriate. Furthermore, the dif- 
fering behaviour away from the extrema implies that for 
consistency we must assume A to be small. One then 
finds dn ~ IQm^viM'^ ~ 0.5GeV^). It is also interesting 
that the effective scale is around M ^ 0.7GeV which is 
well below m„ , and should be cause for concern regarding 
the convergence of the OPE. Nonetheless, one sees that 
the corrections associated with the leading higher dimen- 
sional operators are still quite small. This low scale is 
also the reason we have ignored leading-log estimates for 
the anomalous dimensions, as their status is unclear at 
this scale. A naive application leads to a small correction 
that we shall subsume into our error estimate. 

Extracting a numerical estimate for rf„ from Fig. 2, and 
determining an approximate error, we find the result [20] 



(1-0 ±0-3) 
(500McV)2 



'em* 



(14) 



for the neutron EDM, for which the dominant contribu- 
tion naturally arises from x- 

Comparison with the result of rcf. [6] indicates rather 
good agreement in magnitude, due essentially to the low 



effective mass scale M ^ 700MeV. The relation between 
the calculation presented in this letter and the chiral log- 
arithm estimate will be discussed in more detail elsewhere 
[18], but we note in passing that the sum rule result is 
parametrically enhanced (0(1/A^) vs. 0(l/iV^)) at large 

Combining our result with the recently improved ex- 
perimental bound on d„ [4] we derive, allowing a generous 
50% uncertainty, the limit on theta: 



|6'| < 3 X 10 



-10 



(15) 



which is quite close to previous bounds. 

In conclusion, we have presented a QCD sum rules 

calculation of the ^-induced neutron EDM. This result is 
explicitly tied to a set of vacuum correlators which are 
non-vanishing only in the absence of a U(l) "Goldstone 
boson" . The use of QCD sum rules in the chirally invari- 
ant channel allowed us to unambiguously extract dn{0), 
and independence of the answer from any particular rep- 
resentation of the theta term (7) was checked explicitly. 
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